Eliashberg Analysis of Temperature Dependent Pairing Mechanism in d-wave 
Superconductors: Application to High Temperature Superconductivity. 
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Results are presented for the temperature and frequency dependence of the real and imaginary 
parts of the diagonal self energy for a d-wave superconductor. An Eliashberg analysis, which has 
been successful in recent fitting of superconductor-insulator-superconductor (SIS) tunnel junction 
conductances for Bi2Sr2CaCu208 (Bi-2212), is extended to finite temperatures. The superconduct- 
ing pairing mechanism is assumed to originate in the spin fluctuations of the copper-oxide planes, 
and is modelled by a function incorporating the spin resonance mode measured at an energy of 
approximately 40 meV for optimally doped Bi-2212. The effect of the temperature dependence of 
the spin resonance mode, measured in inelastic neutron scattering (INS), on the finite temperature 
self energies is investigated. 

PACS numbers: 74.20.-z 74.20.Mn 74.25.Jb 74.72.-h 
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Recently, self consistent Eliashberg theory has 
been quite successful in fitting zero temperature 
supcrconductor-insulator-superconductor (SIS) point 
contact tunnel (PCT) junction characteristics in the 
high temperature superconducting cuprate Bi-2212. 1 

The success of that work indicates that the supercon- 
ducting state of the cuprates may be understood to orig- 
inate from a low frequency spectrum of spin fluctuations. 
This interpetation is sometimes referred to as the pairing 
glue scenario 2 , and is most likely to be appropriate for the 
optimal to overdoped regime of the cuprates. For an ex- 
planation of the origin of high temperature superconduc- 
tivity in the underdoped regime, and in proximity to the 
insulating antiferromagnetic phase of the cuprates, the 
RVB singlet picture 3,4 , involving the exchange constant 
J in the t-J model Hamiltonian, provides an explanation 
for the emergence of the high temperature superconduct- 
ing phase. 

A strong coupling analysis has also been used to in- 
terpet both angle resolved photoemission spectroscopy 
(ARPES) 5-7 , and optical conductivity 8,9 . The former 
yields results for the single quasiparticle self energy 
E(w, T). In the latter case, a memory function analysis of 
the optical conductivity results in a frequency dependent 
optical self energy E opt (u;, T) 10 . 

The origin of structure in the electronic spectra of 
the cuprates, measured in such experiments as ARPES 
and optical conductivity, is a topic of much debate. In 
some analyses of ARPES 11 , dispersion anomalies ob- 
served along the nodal line in the copper-oxide plane 
Brillouin Zone, have been interpeted in terms of quasi- 
particle scattering involving phonons, rather than spin 
fluctuations. The energy scales of important features in 
the spin fluctuation spectrum, namely the so-called spin 
resonance mode (at approximately 40 meV), and certain 
phonon modes in optimally doped Bi-2212 make them 
potential candidates for understanding structure seen in 
ARPES. Higher energy spin fluctuations at 80meV have 
been proposed for the observed dispersion anomaly in 
YBaCuO using an Eliashberg analysis 12 . As will be dis- 
cussed later, the absence of a temperature dependence 
in the measured low frequency slope of the real part of 



£(w,T) (Re E(w,T)) in the case of nodal ARPES mea- 
surements on Bi-2212, suggests that the spin resonance 
mode is not involved in nodal quasiparticle interactions in 
Bi-2212, similar to the situation proposed for YBaCuO 12 . 

One motivation for applying a strong coupling, Eliash- 
berg, approach to the cuprates is the success of the 
same approach in conventional superconductors. In con- 
ventional superconductors, such as Pb, a self consistent 
Eliashberg theory, based on the use of a single self energy 
diagram, with no vertex corrections, can be justified by 
Migdal's theorem since the characteristic energy scale of 
phonons relative to the Fermi energy is very small 13 . 

Another important point is that in the case of con- 
ventional superconductors, there exists a detailed, and 
quantitative description of the phonon spectral function, 
denoted by a 2 F(cj), describing the superconducting pair- 
ing mechanism. The analog of a 2 F(ui) for the spin fluc- 
tuation spectrum in the cuprates, denoted by 7 2 x(o;) in 
some work, is not as well understood. The spectral func- 
tion, describing the pairing glue due to spin fluctuations 
is modelled phenomenologically in the present work by 
a prominent peak at the energy of the spin resonance 
mode (ri_R es ) accompanied by a broad background con- 
tinuum extending up to higher energies. 14 The value of 
Jl^es has been measured in inelastic neutron scattering 
(INS) experiments 15 , and is about 40meV for optimally 
doped Bi-2212, which is the focus of the present work. 
The ubiquitous dip feature seen in spectroscopies such 
as tunneling 16 and ARPES 17 is taken as evidence of the 
key role of the spin resonance mode peak in the super- 
conducting state of the cuprates. 

Furthermore, support for the partly phenomenological 
approach of the present work comes from extensive, ab 
initio, Quantum Monte Carlo numerical studies 2, 18-20 on 
the Hubbard and t-J models which indicate a key role is 
played by a low frequency spin fluctuation spectrum in 
determining the properties of the superconducting state 
of these models. 

One of the most unusual aspects of the spin fluctua- 
tion physics in the cuprates is the decrease in intensity 
with increasing temperature of the spin resonance peak at 
tines, along with its disappearence at the superconduct- 
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ing transition temperature Tc- This was first directly 
measured in INS measurements on Bi-2212. 15 Similar be- 
havior has also been extracted from analysis of optical 
conductivity, although in that case the mode appears to 
survive, albeit significantly diminished in intensity, above 
Tc 9 - The emergence of such a peak in the spin fluctua- 
tion spectrum at the superconducting transition may be 
the result of a feedback effect into the electronic spin sus- 
ceptibility due to the onset of the superconducting gap 21 . 

In the pairing glue scenario, the disappearance of 
the resonance mode peak would imply that the super- 
conducting pairing mechanism weakens with increasing 
temperature, very unlike conventional, phonon medi- 
ated, superconductors. The effect of a temperature 
dependent pairing glue on the temperature dependence 
of important quantities such as the superconducting 
gap, and the value of the transition temperature Tc , is 

I 



and 



In the present work 

A(p,u) = A(w,r)cos(20) (3) 

<j) denotes the angular position on the Fermi surface, with 
the cos{2(j>) dependence incorporating a d-wave symme- 
try superconducting state. 

In place of the conventional phonon spectral function 
a 2 F(uj) in the Eliashberg equation formalism, the spin 
fluctuation pairing glue is described by 22 

[c z + c A cos(2(0 - <j>'))]F(w,T) (4) 

The spin fluctuation spectral function F(ui,T) will be 
discussed in detail in the next section. 



one of the subjects of the present work. It also results in 
significant temperature dependent effects in the real and 
imaginary parts of E(w,T) such as a change in slope of 
the ReS(w,T) at low frequencies. Comparison between 
finite temperature model calculations and measurements 
of the temperature dependence of the self energy from 
ARPES and optical conductivity, may serve as an 
important test of the validity of the spin fluctuation 
pairing glue hypothesis for the cuprates. 

Theoretical Formalism 

In the present work, a modified version of the standard 
finite temperature (T ^ 0) Eliashberg equations 13 are 
solved self consistently. The zero temperature (T = 0) 
case of these equations was used in 1 ' 16 . The coupled 
equations to be solved are given by 



(1) 



(2) 

I 

Results 

The results presented here are for optimally doped 
Bi-2212, and are an extension of our successful fitting 
of the SIS tunneling conductance for this material at 
T = 0. 1 The superconducting transition temperature 
is chosen to be 90K in the present work. The values 
of cz = 0.2 and ca = 0.83 were used in reference 
(1) to fit the optimally doped SIS PCT experimental 
conductance. The coupled Eliashberg equations for 
Z(ui,T) and A(uj,T) are solved self consistently Results 
for the diagonal self energy E(w,T) are presented and 
discussed later, where £(w,T) is defined by 

£(w,T) = (1 - Z(u,T))u (5) 

The ReE(w,T) is usually called the mass enhancement 
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factor, and -2 Im£(w,T) is the quasiparticle scattering 
rate (1/t). These two quantities can be measured in 
APRES 5 ' 6 . Optical conductivity experiments 8 yield a 
related optical self energy function as well. The tem- 
perature dependent superconducting gap, usually de- 
noted by A(T), is defined as the frequency uj at which 
u = ReA(w,T). 

The spectral function F(uj,T) describing the pairing 
glue in this work is shown in Figure 1. The intensity of 
the main peak, at the resonance mode frequency 0^ es , 
decreases as the temperature is increased, eventually van- 
ishing at T c . The F(uj,T) curves shown in Figure 1 in- 
corporate the temperature dependence of the resonance 
mode, seen in INS and optical conductivity experiments. 
The resonance mode rides on a weakly temperature de- 
pendent background which extends to higher frequencies. 

Figure 1 is generated using a model for the spectral 
function of the form 

F(u, T) = H(T)F(u)) Res + BG{T)F{uj) bg (6) 

The resonance mode is modelled by a Gaussian F{uj) Res 
with peak magnitude given by H(T) = Hq(1 — 
T/90K) 1 / 2 , where H is the initial peak height at T=0K. 
The temperature dependence of the peak height is similar 
to the BCS mean field A(T) behavior, and also consis- 
tent with the temperature dependent mode height seen 
in INS experiments. 15 

The background continuum is modeled using a broad 
temperature independent Gaussian function of frequency 
oj denoted by F{uj) bg multiplied by a temperature de- 
pendent factor BG(T) = 1 + aT/90K, where a is ad- 
justed to yield A(T) = at T = 90K in the self consis- 
tent solution of the Eliashberg equations. 



with INS measurements is not understood. However, if 
the temperature dependence seen in optical conductivity 
was used in place of the complete mode disappearance as 
used in the current work, the behavior discussed below in 
the 2A(0)/fc B T c and the self energy E(w,T) would still 
be present, albeit quantitatively different. 

In one analysis (see Fig (18) of reference (10)) of 
ARPES experimental data, the resonance mode peak 
sits on a relatively large, temperature independent, back- 
ground continuum. In that case, the decrease in the peak 
height, as the temperature is increased up to Tc results 
in a correspondingly smaller reduction in the pairing glue 
coupling strength compared to the current results pre- 
sented here. 

Finally, the issue of whether the area under the F(u,T) 
curve should be conserved, as the temperature changes, is 
unclear. This would require that area lost from the main 
peak would be redistributed to higher frequencies. How- 
ever, even if this requirement were imposed, there would 
still likely be a reduction in the strength of the pairing 
glue (see equation (7) ahead) with increasing tempera- 
ture. 

The consequences of a complete disappearance of the 
resonance mode peak are investigated in this work, but 
uncertainties in the frequency dependence of the pairing 
glue spectral function should be kept in mind. 

One of the interesting features of Bi-2212 (optimally 
doped) is the high value of the ratio 2 A(0) / ksTc , which 
is observed to be 9.3. Table 1 lists values from reference 
(23) of measured Tc and zero temperature gap values 
(A(0)) and the corresponding 2A(0) /k B T c . 

TABLE I: T c and A(0) values from reference (23) 



T C (K) A(0) (meV) 2A(0)/fc B T c 




E (meV) 



FIG. 1: Evolution of pairing glue spectral function F(uj,T) in 
equations (1) and (2) for various temperatures up to T c . 

In the present work, the resonance peak at £ln es disap- 
pears completely at Tc — 90K. 15 This is consistent with 
INS measurements of the temperature dependence of the 
spin resonance mode intensity. In a recent analysis of 
optical conductivity on Bi-2212 9 , the height of the reso- 
nance mode decreases by a factor of two at the supercon- 
ducting transition temperature, with a smaller peak sur- 
viving into the normal state. The origin of this difference 
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The disappearance of the resonance mode at the mea- 
sured Tc may be an important contribution to the large 
2A(0)/fc B T c ratios in Table 1. 

To illustrate the effect of assuming a temperature inde- 
pendent pairing glue, Figure 2 shows the resulting A(T) 
from using the T = OK F(lu, T) from Figure 1 for all 
finite temperatures in the Eliashberg equations. The 
calculated T c for this case is 153K as seen in Figure 
2, much higher than the measured 90K value for opti- 
mally doped Bi-2212. The value of the ratio 2A(0) /k B T c 
is 5.7, much less than the observed value of 9.3. This 
same discrepancy will occur in the case of the other Bi- 
2212 doping levels previously analysed 1 . This observa- 
tion implies that a temperature dependent weakening of 
the strength of the pairing glue, and a consequent re- 
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duction in the transition temperature Tc, may be at the 
origin of the large ratio values shown in Table 1 Fur- 
thermore, the trend seen in Table 1 for the 2A(0)/A;bT c 
ratio to decrease with overdoping is also consistent with 
the decreasing prominence of the resonance mode in the 
pairing glue spectrum with overdoping 1 . 




20 40 60 80 100 120 140 160 
T(K) 

FIG. 2: A(T)/A(0) vs T using the F(u),T) for T = from 
Figure (1). T c is 153K. 

Figures (3) and (4) show the Re£(w, T) and ImS(w, T) 
for the case where F(u>,T) is kept at its T = shape 
for all temperatures up to the transition temperature 
of 153K. As temperature increases, the peak structures 
in the self energy functions are smeared out due to the 
Fermi and Bose functions in the Eliashberg equations. 
It is worth noting that, in the d-wave superconduct- 
ing state, the main peak in the Re£(w,T = 0) is lo- 
cated at approximately 75 meV, slightly below the value 
A(0) + n Res = 82 meV (using A(0) = 38meV and 
Flues = 44 meV 1 ), which would be expected for an s- 
wave superconductor. 

A measure of the pairing glue coupling constant in- 
volved in determining £(u;,T), which is denoted by A in 
the present work, can be extracted in two ways. One is 
by analogy with the electron-phonon coupling constant 
in conventional superconductors, by use of 

oo 

A = 2c z J F(u,T)/u (7) 
o 

A second measure of A is from the slope, at low uj, of 
ReE(w,T). Results for A from both of these procedures 
are tabulated in Table II for the case of a temperature 
dependent pairing glue, which will be discussed later. For 
the current situation where the F{uj,T) is assumed to 
stay constant at its T = OK shape, the slope of ReS(a;, T) 
in Figure (3) is constant at low u>, and consequently A 
does not change with temperature. 

It is interesting to compare the temperature evolution 
seen in Figures (3) and (4) with that of Pb for which the 
Re£(u;,T) and lmS(u,T) are shown in Figures (5) and 
(6). The lead gap A(0) = 1.35meV is almost a factor of 3 
less than the first phonon peak of approximately 4.4meV 
in the lead phonon a 2 F(uj) spectral weight, which itself 
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FIG. 3: Temperature dependence of the real part of the self 
energy S(w) using the T — F(lj,T) from Figure (1) for all 
temperatures shown in this figure. 
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FIG. 4: Temperature dependence of -ImE(w) for same con- 
ditions as in Figure (3). 



is temperature independent. The result is a relatively 
undramatic temperature evolution in the Pb ReS(cj,T) 
and ImE(w,T) compared with that seen in Figures (3) 
and (4) where the value of A(0) and flR es are comparable. 

Figures (7), (8) and (9) show results for A(T), 
Re£(w,T) and ImE(w,T) for the case of the tempera- 
ture dependent F{lo,T) shown in Figure (1). The su- 
perconducting gap A(T) now goes to zero at 90K, due 
to the choice of a = 0.6 in the background continuum 
BG(T) function. The T = value for A(0) = 38meV 
is in agreement with measured value from SIS PCT 
measurements 23 , and so the correct 2A(0)/&bTc; of 9.3 
is reproduced. The overall temperature dependence of 
A(T) in Figure (4) resembles that of the resonance mode 
peak height temperature dependence used for H(T). 

Figures (8) and (9) show the corresponding temper- 
ature dependent real and imaginary parts of E(w,T). 
There are significant differences between these results 
and those shown in Figures (3) and (4). The ReE(w, T) of 
Figure (8) show a much more rapid decrease in the main 
peak at approximately Qn es + A(T) which reflects the 
decrease in the intensity of the resonance mode with in- 
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FIG. 5: Temperature evolution of ReE(aj,T) for lead (Pb) 
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FIG. 6: Temperature evolution of -Im£(w,T) for lead (Pb) 



creasing temperature. The dip feature also fills in more 
rapidly compared to the temperature evolution of Fig- 
ure (3). Furthermore, there is now a noticeable decrease 
in the slope of Re£(u;,T) at low frequencies, reflecting 
the decrease in the effective coupling constant (A) as the 
resonance mode disappears. This is illustrated in the ac- 
companying Table II. 

In Figure (9), the disappearance of the spin resonance 
mode peak with increasing temperature results in a sig- 
nificant reduction in Im£(u;,T) in the energy range 50 
meV to lOOmeV. This energy range brackets the energy 
fines + A(T) where the mode shows up in such quanti- 
ties as the self energy in Eliashberg theory. A significant 
reduction in the measured width of the MDC peak in 
the same frequency range, measured by ARPES, could 
search for this effect. 

A limited amount of experimental measurements of the 
real and imaginary single particle H(u>,T), have been ex- 
tracted from ARPES MDC curves on Bi-2212 5 ' 6 . Optical 
conductivity experiments on Bi-2212 8 have also yielded 
an optical self energy £ opt (w,T), along with a tem- 
perature dependent spectral function for the underlying 
modes interacting with the quasiparticles. These anal- 
yses appear to reveal several different contributions to 
the extracted self energies. Separating out different con- 




20 40 60 80 100 

T(K) 

FIG. 7: Temperature dependent superconducting gap A(T) 
for optimal doped Bi-2212 for the temperature dependent 
pairing glue spectra F(ui,T) shown in Figure (1) 



TABLE II: Coupling strength of optimally doped Bi2212 for 
various temperatures up to T c (from Figure (8)) cz = 0.2 (see 
text) 



T(K) \ 



0.01 0.46 0.59 

36 0.45 0.57 

54 0.44 0.54 

72 0.41 0.49 

85.5 0.36 0.41 

90 0.23 0.29 



tributions to the real and imaginary components of the 
quasiparticle self energy has been the subject of recent 
work (see Figure (10) of reference (10)) and reference 
(24). 

In reference(8), a contribution is present in 1/t = - 
2ImE opt (aj, T) which increases linearly with energy, and 
is temperature dependent. The linear behavior with 
energy in 1/t has been interpeted as evidence for the 
Marginal Fermi Liquid hypothesis 25 . Furthermore, a 
large energy and temperature dependent background is 
measured in the Re£(w, T) in reference (6) which is sub- 
tracted off to yield the temperature dependent contri- 
bution to ReE(u;,T). These contributions complicate 
interpetation of the ARPES and optical conductivity 
data, and their comparison with the present model. In 
the analyses of ARPES and optical conductivity experi- 
ments, a background continuum extending up to 400meV 
is present in the I 2 x{u) curves (see Figures (18) and (31) 
of reference (10)), in contrast to the pairing glue spectral 
curves shown in Figure (1) of this work where the high 
energy continuum extends up to 175 meV approximately. 
The linear increase in 1/t with increasing energy in ref- 
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E (meV) 

FIG. 8: Temperature dependence of ReE(oj) for the temper- 
ature dependent F(u>,T) spectra shown in Figure (1). 



T=0.01K 

T=36.0K 




E (meV) 

FIG. 9: ImE(w) for same conditions as in Figure (8). 



erence (8) most likely results from this. 



Bearing these issues in mind, figures 3(b) and 4(b) of 
reference (6), which show the MDC linewidth and the 
Re£(w, T) for nodal line ARPES of Bi-2212 respectively, 
for different temperatures, do not appear to display the 
rapid decrease in linewidth or a discernible change in 
slope for Re£(u;,T) at low frequencies, that are shown 
in Figures (9) and (8) of the present model calculations. 
This could be interpeted as evidence that the nodal Bi- 
2212 quasiparticles are not coupling to the 40meV tem- 
perature dependent spin resonance mode. The zero su- 
perconducting gap for these quasiparticles would point 
to a central role for the 40meV mode in the supercon- 
ducting pairing in Bi-2212. Experimental studies of the 
anti-nodal self-energies would be of interest in following 
up on this issue. 



In conclusion, results have been presented for an finite 
temperature Eliashberg calculation of a d-wave supercon- 
ducting case with temperature dependent pairing glue. 
The calculated real and imaginary diagonal self energies 
may prove a useful probe of the pairing glue description of 
high temperature superconductors when compared with 
ARPES and optical conductivity experiments. 



The authors wish to acknowledge useful conversations 
with J. F. Zasadzinski 



1 O. Ahmadi, L. Coffey, J. F. Zasadzinski, N. Miyakawa and 
L. Ozyuzer et al, Phys. Rev. Lett. 106, 167005 (2011) 

2 T.A. Maier, D. Poilblanc, and D.J. Scalapino, Phys. Rev. 
Lett., 100, 237001 (2008) 

3 P. W.Anderson, Science,316, 1705 (2007) 

4 T. M. Rice, K. Y. Yang, and F.C. Zhang, Rep. Prog. Phys. 
75, 016502 (2012) 

5 P. D. Johnson, T. Valla, A.V. Fedorov, Z. Yusof, A.R. 
Moodenbaugh, G.D. Gu, N. Koshizuka, C. Kendziora, Sha 
Jian, and D.G. Hinks, Phys. Rev. Lett. 87 177007 (2001). 

6 W. Zhang, G. Liu, L. Zhao, H. Liu, J. Meng, X. Dong, W. 
Lu, J. S. Wen, Z. J. Xu, G. D. Gu, T. Sasagawa, G. Wang, 
Y. Zhu, H. Zhang, Y. Zhou, X. Wang, Z. Zhao, C. Chen, Z. 
Xu, and X. J. Zhou, Phys. Rev. Lett. 100, 107002 (2008). 

7 T. Valla, T. E. Kidd, W. G. Yin, G. D. Gu, P. D. Johnson, 
Z. H. Pan, and A.V. Fedorov, Phys. Rev. Lett. 98, 167003 
(2007) 

8 J. Hwang, T. Timusk, and G. D. Hu, Nature (London), 
427, 714, (2004) 

9 J. Hwang, T. Timusk, E. Schachinger and J. P. Carbotte, 
Phys. Rev. B 75 144508 (2007) 

For a recent review and references therein see J.P.Carbotte, 
T. Timusk, J. Hwang, Rep. Progress Physics, 74, 066501 
(2011) 



11 S. Johnson, I. M. Vishik, W.S. Lee, F. Schmitt, S. Uchida, 
K. Fujita, S. Ishida, N. Nagaosa, and Z.X. Shen, arxiv.org 
cond-mat 1101.1302. 

12 T. Dahm, V. Hinkov, S.V. Borisenko, A. A. Kordyuk.V.B. 
Zabolotny, J. Fink, B. Buchner, D.J. Scalapino, W. Hanke, 
B. Keimer, Nature, 5, 217 (2009) 

13 D. J. Scalapino, J.R. Schrieffer, and JW. Wilkins, Phys. 
Rev. 148, 263 (1966) 

14 E. van Heumen, E. Muhlethaler, A. B. Kuzmenko, H. 
Eisaki, W. Meevasana, M. Greven, and D. van der Marel, 
Phys. Rev. B79 184512 (2009) 

16 H.He, Y. Sidis, P. Bourges, G.D. Gu, A. Ivanov, N. 
Koshizuka, B. Liang, C. T. Lin, L.P. Regnault, E. Schoen- 
herr, and B. Keimer, Phys. Rev. Lett. 86 1610 (2001) 

16 J. F. Zasadzinski, L. Ozyuzer, L. Coffey, K. E. Gray, D. 
G. Hinks, and C. Kendziora Phys. Rev. Lett. 96, 017004 
(2006) 

17 Z.-X. Shen, D. S. Dessau, B. O. Wells, D. M. King, W. E. 
Spicer, A. J. Arko, D. Marshall, L. W. Lombardo, A. Ka- 
pitulnik, P. Dickinson, S. Doniach, J. DiCarlo, T. Loeser, 
and C. H. Park, Phys. Rev. Lett. 70, 1553 (1993 

18 B. Kyung, D. Senechal, A. M.S. Tremblay, Phys. Rev. B 
80, 205109 (2009). 

19 T. Maier, A. Macridin, M. Jarrell, and D. J. Scalapino, 



7 



Phys. Rev. B 76 144516 (2007). 

E. Khatami, A. Macridin, and M. Jarrell, Phys. Rev. B 
80, 172505 (2009) 

T. Dahm, D. Manske, and L. Tewordt, Phys. Rev. B. 58, 
12454 (1998) 

K.A. Musaelin, J. Betouras, A.V. Chubukov, and R. Joynt, 
Phys. Rev. B 53 3598 (1996)). 

J. F. Zasadzinski, L. Ozyuzer, N. Miyakawa, K. E. Gray, 



D.G. Hinks, and C. Kendziora, Phys. Rev. Lett. 87, 067005 
(2001) 

S. D. Conte, C. Giannetti, G. Coslovich, F. Cilento, D. 
Bassini, T. Abebaw, F. Banfi, G. Ferrini, H. Eisaki, M. 
Greven, A. Damascelli, D. Van der Marel, F. Parmigiani, 
Science, 335 1600 (2012) 

C. M. Varma et al, Phys. Rev. Lett. 63, 1996 (1989) 



